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ABSTRACT 

We present a model for the cosmological evolution of quasars (QSOs) under the 
assumption that they are powered by massive accreting black holes. Accretion flows 
around massive black holes make a transition from high radiative efficiency (~ 10%) 
to low efficiency, advection-dominated flows when M/M^dd falls below the critical rate 
~ 0.3a 2 ~ 10 -2 where M is the mass accretion rate, M^dd oc M is the usual Eddington 
rate with the nominal 10% efficiency, and a(< 1) is the dimensionless viscosity 
parameter. We identify this transition with the observed break at a redshift ~ 2 in 
the QSOs' X-ray luminosity evolution. Growth of black holes through accretion could 
naturally lead to such a transition at a critical redshift z c ~ 1 — 3, provided that most 
of high redshift QSOs appear with near Eddington luminosities at z ~ 3 — 4 and the 
accretion rates decline over the Hubble time in a roughly synchronous manner. Before 
the transition, the QSOs' luminosities (with a high efficiency) slowly decrease and after 
the transition at z c , the QSO luminosities evolve approximately as oc (1 + z) K ^ where 
K{z) gradually varies from z = z c to z ~ around K ~ 3. The results depend on the 
details of the QSO X-ray emission mechanism. We discuss some further implications. 

Subject headings: accretion, accretion disks — cosmology: theory — quasars: general 



Introduction 



The cosmological origin and evolution of quasars (QSOs) has remained an outstanding issue 
in astrophysics. QSOs first appear at a high redshift z > 4 and their evolution appears very slow 
between z ~ 3 and the critical redshift z c ~ 2. Below z c , bright QSOs seem rapidly disappearing 
(for a review, see Hartwick & Shade 1991). In the optical, the evolution is often described as 
pure luminosity evolution (Mathez 1976,1978, Weedman 1986) with the luminosity decreasing as 
(1 + z) K with K ~ 3 — 4 depending on the cosmological model (Marshall 1985, Weedman 1986, 
Boyle et al. 1988, Hartwick & Shade 1991). In the X-ray, a similar evolution is seen although the 
optical evolution appears somewhat steeper (Maccacaro et al. 1991, Delia Ceca et al. 1992, Boyle 
et al. 1993, 1994, Page et al. 1996). 

In the conventional picture (for a review, see e.g. Rees 1990, Frank et al. 1992), QSOs are 
powered by massive accreting black holes which grow through accretion of gas. It remains largely 
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unexplained why QSOs are suddenly quenched at z < 2 (cf. Turner 1991, Fukugita h Turner 
1996). A tempting possibility that low z Seyfert galaxies may be low luminosity remnants of QSOs 
is ruled out simply due to the fact that the former are more abundant than the latter at least by 
an order of magnitude (Huchra & Burg 1992). Why do QSOs suddenly turn off? What drives the 
evolution of the QSO luminosity function? What are the present remnants of bright QSOs? We 
address these issues. 

We assume that engines of QSOs are indeed massive accreting black holes. The accretion 
flows have often been modeled by a geometrically thin, optically thick disk (e.g. Frank et al. 1992, 
Pringle 1981) and the so-called blue bumps have been interpreted as emission from such disks (e.g. 
Sun & Malkan 1988, Wandel & Petrosian 1988, Laor & Netzer 1989, and references therein). For 
the optical QSO evolution, some phenomenological prescriptions have been proposed (Heisler & 
Ostriker 1988, Koo & Kron 1988, Boyle et al. 1988,1993, Pei 1995). For the X-ray evolution, it has 
been shown that the pure luminosity evolution can account for the observed evolution (Boyle et 
al. 1993,1994, Page et al. 1996). Although there have been various models on the possible driving 
mechanism for the QSO cosmological evolution (e.g. Caditz et al. 1991, Small & Blandford 1992, 
Haehnelt & Rees 1993), a definite conclusion has not been reached. 

For a typical QSO luminosity L ~ 10 46 erg/s, a black hole of mass m = M/10 8 Mq accreting 
at the Eddington rate with the nominal radiative efficiency i] n = 0.1, 

M EM = L Edd /{r] n c 2 ) = AttG M / r] n n es c = [1.4 x I0 26 g/s]m, (1-1) 

where K es is the electron scattering opacity, would exponentially grow on a time scale 

-l 



tEdd = 



M Edd /M = [4.5 x IOVKWO.1). (1-2) 



This suggests that even when accretion-powered QSO activity is short-lived (< 10 9 yr), massive 
black holes residing in most luminous QSOs (L > 10 46 erg/s) could grow to M > 1O 9 M . We 
suggest that this rapid growth itself may directly affect the QSOs' evolution. The basic idea is that 
the massive black holes, accreting at a scaled accretion rate m = M /Mem below a certain critical 
rate m c ~ 10~ 2 (Narayan & Yi 1995b), emit radiation with an efficiency i] cx m <C r\ n . Although 
the idea of the low accretion efficiency (Rees et al. 1982) has occasionally been mentioned, the lack 
of a self-consistent model has been a major obstacle. Since Mem °^ M increases due to accretion 
(even when M remains nearly constant), m cx M/M decreases and the QSOs' luminosity decrease 
could be a direct result of the rapid mass accretion. We suggest that the transition in accretion 
flow from high efficiency (rj ~ r] n = 0.1) to low efficiency (rj <C r] n ) flows could naturally explain 
the major properties of the X-ray luminosity evolution of QSO population (Yi 1996). 



2. Cosmological Evolution of QSO Luminosity: Transition in Accretion Flows 
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The idea of transition in accretion flow hinges on the recently proven stability of the 
advection-dominated flows (Narayan & Yi 1994, 1995ab, Abramowitz et al. 1995, Chen et al. 
1995). As the accretion rate falls, the density of the accreting plasma drops to a point where 
the Coulomb ion-electron energy transfer becomes inefficient compared with the rate of viscous 
heating on ions (Rees et al. 1982). In the absence of other couplings, the accretion flow becomes 
a two temperature plasma with a very low radiative cooling efficiency (Narayan & Yi 1995b, 
Rees 1982) and most of accretion energy is lost into black holes through advection. Optically 
thin (magnetized) advection-dominated flows emit synchrotron (in radio) and Comptonized 
synchrotron and bremsstrahlung emission (in X-ray) resulting in characteristically flat spectra. 
Detailed spectral calculations and results have been successfully applied to low luminosity systems 
such as Galactic black hole systems (Narayan et al. 1996), the Galactic center (Narayan et al. 
1995), and the extragalactic system NGC 4258 (Lasota et al. 1996). In the case of the Galactic 
center, the estimated accretion efficiency rj ~ 10 _5 a _1 and in NGC 4258 7] ~ 10 _4 a _1 where 
a ~ 10 -2 — 1 is the dimensionless viscosity parameter (e.g. Pringle 1981, Prank et al. 1992). 
Fabian and Rees (1995) have also suggested that dim galactic nuclei in large elliptical galaxies 
(due to low efficiency advection-dominated flows) contain remnants of luminous QSOs. 

An accretion flow enters the low efficiency regime at a critical accretion rate, m c = M c j 'M^dd-, 
which is primarily determined by a{< 1) (Narayan & Yi 1995b), 

m c « M c /M Edd = 0.3a 2 (2-1) 

for R < 10 3 R S where R is the radial distance and R s = 2GM/c 2 . For R > W 3 R S , 

m c ^0.3a 2 {R/l0 3 R s )- 1/2 . (2-2) 

Since most of the accretion energy is dissipated at R < 10 3 R S , the critical rate given in eq. (2-1) is 
relevant. When m > m c , the accretion flow takes the form of high efficiency (i.e. r] ~ r] n = 0.1) or 

L = 7] n Mc 2 . (2-3) 

For m < rh c , the low efficiency is approximately described by the luminosity 

L « 30m x L Edd (2-4) 

with x ~ 2 — 2.2. We take a = 0.3 or m c ~ 0.03 and x = 2.2 for our discussions (Narayan & Yi 
1995b). Unless the accretion rate rapidly increases toward the present epoch, the growth of black 
holes always tend to decrease m = M / MEdd oc M/M, which naturally drives the transition. The 
overall evolution takes the form of the the pure luminosity evolution. 

The direct comparison with the observed X-ray evolution is complicated by the uncertainties 
in X-ray emission from luminous QSOs. The conventional thin disks (with high efficiency) have 
the effective disk temperatures < 10 6 K and hence cannot emit X-rays efficiently (e.g. Frank 
et al. 1992). We assume that a significant fraction of the bolometric luminosity is emitted as 
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X-rays (presumably from accretion disk coronae) and that rj n = 0.1 represents this uncertainty. 
The advection-dominated hot flows differ from the thin disks as the former emit hard X-rays up 
to several 100 keV. As long as most of the accretion flow is advection-dominated, the emission 
spectra are characteristically flat, i.e. luminosity per decade of frequency is roughly constant (e.g. 
Narayan et al. 1995, Lasota et al. 1996). The X-ray luminosity can be approximated as a fixed 
fraction of the total luminosity. In our discussions, we simply assume that the observed X-ray 
luminosity is proportional to the bolometric luminosity. We caution that our assumption could 
introduce serious errors if the X-ray to bolometric luminosity ratio is a complex function of some 
yet unknown parameters such as L, M, or rh. For instance, if the accretion flow contains a thin 
disk in the outer region (eq. (2-2)), the resulting emission spectra could be significantly affected. 
In this sense, it is uncertain how our suggestion on the X-ray luminosity evolution could be applied 
to the optical evolution (see below). 

We begin with a simplest scenario; (i) The black holes are initially accreting at their 
Eddington rates (cf. Padovani 1989). (ii) The accretion rates are sufficiently slowly decreasing 
from the initial epoch at cosmic time t = U (or z = Zi). For a flat universe (Q.q = 1) with no 
cosmological constant (ft a = 0), the cosmic time t at z is given by 



H (t — ti 



(l + z 



1-3/2 _ 



(i+^r 3/2 



(2-5) 



where U = t(z = z^). When M « constant, or 



M ss Mi + M(t - ti 



(2-6) 



where Mj = M(z = z{), 



rh oc M 1 oc 



i +3!fc( (i+ ^ ,/2 - (i+z - ) - 3/2 ) 



(2-7) 



where H Q = 50km/s/Mpc is the assumed Hubble constant. While rh > rh c , the observed slow 
luminosity evolution in X-rays and optical could be purely due to the gradual decrease of M (eq. 
(2-3)). When rh < rh c , the luminosity evolves as (eq. (2-4)) 



L oc 



-x+l 



(2-8) 



In the limit that the accreted mass over Hubble time is larger than the initial mass, i.e. 
Mj H » M h 



L oc (1 + zp- 1 )/ 2 (1 + Zi fl* - (1 + zfl 2 



-x+l 



(2-9) 



Similarly, for an empty universe (Qq = Cl\ = 0), 



H (t-U) = [(1 + z)- 1 ' 2 - (1 + Zi)- 1 '* 



(2-10) 
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and 



L oc 



1+ A^(< 1 + ^ 1/2 -< 1+ ^ 1/2 ) 



M ( if„ 



-X+l 



(l + z 4 ) 1/2 -(l + ^) 1/2 



-x+l 



or in the limit M/H D > M, 

Lo,{l+z)^' 2 

The simple behavior of luminosity as a function of redshift is often approximated as 

L(z) oc (l + z) K , 



(2-11) 

(2-12) 
(2-13) 



with a constant K. In our scenario, -ftT(z) is a function of 2 & Zj and depends on the cosmological 
model. We get for the flat universe, 



K(z) 



dlnL(z) 
dln(l + z) 



3fs - 1) 



1 + 



1 



and for the empty universe, 



K(z) 



x 



1 



1 + 



((1 + ^/(1+^)3/2-1 



1 



(2-14) 



(2-15) 



((l + ^)/(l +2 ))l/2_l_ 

where eqs. (2-9), (2-12) have been used. The redshift dependence of K is shown in Fig. 1. In the 
two cosmological models considered, which cover various other possible models, the evolution at 
z < zi is shown to give a wide range of K which includes K ~ 2 — 3. K varies from K > 3 to 
K < 3 as z decreases away from Zj. The derived expressions for K{z) become invalid for z — > Zj as 
the accretion flows become those with a high efficiency in our scenario. The derived values of K 
are interestingly close to those obtained by Boyle et al. (1994) and Page et al. (1996) using the 
recent X-ray observations (see below). 

Another important aspect of the QSO evolution is the sudden cutoff in the luminosity 
evolution near z = z c ~ 2 which we identify as the critical redshift at which m = m c first occurs. 
For the flat universe with no cosmological constant, we get 

-1 



n o 2 r] n K es c M 
m c = 0.3a = — — — — — 
4vrG Mi 



1 + 



2M 



( l + Zc) -3/2_ (1 + z . r 3/2 



If the initial accretion rate is a fraction <5(< 1) of the initial Eddington rate or 

M/M t = 5r E \ d = (4.5 x lOVrWO-l) -1 *, 

-2/3 



(2-16) 



(2-17) 
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/ tEdd \ 
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(2-18) 



where t 



age 



2tff/3 = 2/3H is the age of the Universe. Similarly, for the empty universe 



(Qo = n A = o) 



l + z c 



tEdd 



''age 



1 

rh c . 



+ (1 + Zi 



-1/2 



(2-19) 
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where t age = tn = 1/H Q . Unless S < rh c ~ 0.03 (which is unlikely given the high luminosities of 
QSOs at high redshift z = Z{ and a relatively short time scale allowed for growth of black holes 
from z > Zj), z c only weakly depends on <5 ~ 1. Assuming a = 0.3, 6 = 1, and H a = 50km/s/Mpc, 
we get z c 1.6(1.8) for Zj = 3(4) in the flat universe and z c 1.9(2.6) for = 3(4) in the empty 
universe. Increasing H to i? = 65km/s/Mpc would change z c roughly by ~ 10 — 20%. For 
rii c <C~ 3 x 10 -2 or a < 0.3, however, z c could be much lower than the observed z c ~ 2 unless 
the initial accretion rate is already significantly sub-Eddington. Our proposed evolutionary model 
suggests a possible explanation for the cutoff at z c ~ 2. A definite final answer requires fine 
details of the cosmological model and the initial conditions. Nevertheless, it remains valid that the 
transition driven by mass growth could always result in rapid evolution of QSOs regardless of any 
specific cosmological scenario such as the present one. 

Now we consider luminosity evolution in general cosmological models relaxing M = constant. 
We adopt the mass accretion rate of the form 

M = Miexp(-(t-ti)/t evol ) (2-20) 

with the evolutionary time scale as t evo i = [3t age ((3 < 1). In order to see the effects of the 
cosmological model, we consider four different cosmological models (e.g. Kolb & Turner 1990); (i) 
an open universe fin = 0.4, S7a = for which t age = 0.779(l/-ff o ) and 



t(z) 



1 



2H Q (1 - f) )3/2 

2(i- n ) 1 / 2 (i + n z) 1 / 2 
n (i + z ) 



fi (l + z) + 2(1 - fi ) + 2(1 - f\i) 1/2 (l + fV) 1/2 



(ii) an empty universe f^o = = (eq. (2-10)), (iii) a flat universe f^o = 1, = (eq. (2-5)), 
(iv) a flat universe Oq = 0.4, S7a = 0.6 for which t age = 0.888(l/i/ o ) and 



t(z) = 



3H (1 - Q ) 1/2 



In 



On 




(2-22) 



In each cosmological model, we need to specify initial conditions for QSO evolution. Our simple 
initial conditions could be partly motivated by the optical observation that the QSO luminosity 
function is already established by z ~ 3 — 4 and changes little from z ~ 3 to z ~ 2 (Hartwick & 
Shade 1990). The present model simply posits that the massive black holes arrive at the peak of 
their activities (z ~ 3 — 4), in a nearly synchronous manner, accreting from surrounding gas at a 
rate close to Mem- 

In Fig. 2, we show evolution of QSOs in the four cosmological models marked by numbers. We 
have taken Z{ = 3.5 in this example. In our model, the luminosity evolution is obviously scale-free 
for our chosen initial conditions. That is, the luminosities of QSOs are simply determined by the 
black hole masses with an identical redshift dependence. The most significant uncertainty is the 
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possible wide scatter in the initial conditions (i.e. Zi, 5, t evo i, etc.). At z > z c , the evolution is slow 
(eqs. (2-3), (2-20)) and the time scale of evolution is set by t evo [ as % = 0.1 is maintained. After 
the transition at z = z c , the difference in the power-law among the shown models is small. It is 
encouraging that the observed constant K ~ 3 between z = and z = z c ~ 2 could be accounted 
for with small effects from the cosmological model. For the same set of parameters and the initial 
conditions, z c significantly increases from the flat universe to the empty universe. Such an effect 
is, however, indistinguishable from that due to uncertainties in rh c or a. t evo i <C t age or <C 0.5 
could result in much faster evolution than the observed. A faster decrease of M near z = would 
inevitably lead to the further acceleration of evolution which might be directly observable as 
deviation from the observed power-law evolution with a constant K. 

X-ray observations have established a well defined evolutionary trend although details differ 
in several analyses. Delia Ceca et al. (1992) and Maccacaro et al. (1991) obtained K = 2.56 
and Boyle et al. (1993) derived K = 2.75 - 2.80 for the empty universe and K = 2.76 - 2.73 
with z c ~ 2 for the flat universe with no cosmological constant. More recently, Boyle et al. 
(1994) analyzed a larger sample and re-derived K = 3.34 [z c = 1.79) for the empty universe and 
K = 3.25 (z c = 1.60) or K = 3.17 (z c = 1.48) for the flat universe. Page et al. (1996) note that 
no evolution for z > 1.8 is entirely compatible with the observed evolution with K = 2.35 — 2.94 
& z c = 1.41 — 1.82 depending on the cosmological model. Despite lacking details of evolution, 
our proposed model is largely consistent with the observations. The X-ray evolution is similar to 
the optical evolution except that the latter is slightly faster than (cf. Boyle et al. 1988, Boyle et 
al. 1994). In the optical, the derived K depends rather sensitively on the assumed spectral slope 
which could be a major source of uncertainty (see below). 

Caditz et al. (1991) has suggested that optical evolution could be explained by the accretion 
disk model where the monochromatic luminosity at a given frequency evolves rapidly as the fixed 
frequency crosses the Wien cutoff frequency due to the increase of the black hole mass through 
accretion. Interestingly, such an explanation also results in a luminosity evolution of the form 
L (x rh 2 (cf. eq. (2-4)) in the optical. In other words, the observed (1 + z) K (K ~ 3 — 4) luminosity 
evolution could be the result of L oc rh 2 arising from the two entirely different origins. In our 
model, the long wavelngth emission is unclear due to the uncertain nature of the outer region of 
the flow where both a thin disk and an advection-dominated flow are allowed in principle. For 
instance, according to eqs. (2-1), (2-2), there could be an outer thin disk while the inner region is 
fully advection-dominated when 

0.3a 2 (R/10 3 R s y 1/2 < rh(R) < 0.3a 2 (2-23) 

at a distance R > 10 3 i? s . If the outer disk continues to exist at z < z c , we can draw a rough 
estimate on the bolometric luminosity evolution at long wavelengths by directly applying eq. (2-2). 
Unless the outer disk extends to a region well inside ~ W 3 R S , however, the optical luminosity is 
expected to be very low and the following evolutionary trend is not directly applicable to brightest 
QSOs. As rh oc M/M decreases, the inner region of the disk first becomes advection-dominated 
while the outer thin disk, emitting in the optical, is maintained. As rh drops further, a larger 
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fraction of the inner disk becomes advection-dominated and the thin disk is pushed outward. The 
inner radius of the outer thin disk is 

0-3« 2 {-^rV 2 = m oc -L (2-24) 



\W 3 R S J M 
or 

R in as [2.7 x 10 13 cm]mm" 2 (2-25) 

for M ~ constant and a = 0.3 as before. The bolometric disk luminosity from the outer disk 
would decrease as 

GMM MM _ 2 . . 

as M increases. 



or 



The disk emission would have the evolving peak temperature 

/ mm\ 

-1/4 / ^ X7/4 



^ ~ ( ^} oc M- 2 oc L disk . (2-27) 



T dls ^[2.6xl0 3 K](^)' (^) . (2-28) 
The peak wavelength is 

The emission from the outer disk is mostly in the infrared regime and the optical emission 
from the disk will occur in the Wien tail of the emission at ~ T<n s k- The monochromatic disk 
luminosity at a fixed optical wavelength X opt < X p would rapidly decline as X p oc M 2 increases and 
X pt A p . In this case, the observed optical-UV emission is dominated by that from the inner 
advection-dominated flow and the optical-UV evolution is similar to the X-ray evolution. At a 
longer infrared wavelength \m > \ p , the monochromatic luminosity 

Lm oc T disk /\j R oc M 4 (2-30) 

where we have used the Rayleigh-Jeans law. In other words, the infrared luminosity at Xjr 3> A p 
would increase as oc M 4 until X[r ~ A p is reached. When A p > X[r occurs subsequently, the 
infrared luminosity also rapidly declines. One of the major uncertainties in this qualitative 
prediction is the effects of dust and irradiation from the central X-ray emission region which could 
seriously affect the long wavelength emission from the outer region. 
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3. Implications 

If the QSO evolution is indeed due to a single population, the entire QSO population shifts 
to a mass range > 10 9 M© by z = 0, which suggest that massive elliptical galaxies are most likely 
hosts for QSO remnants (Fabian & Canizares 1988, Fabian & Rees 1995). Such massive black 
holes have not been clearly identified with an exception of M87 (Kormendy and Richstone 1996). 
When the radiative efficiency is a sensitive function of the mass accretion rate, the mass estimates 
based on r] n ~ 0.1 could be seriously misleading (Narayan & Yi 1995b, Fabian & Rees 1995). In 
our model, the mass of a black hole at z = 0, based on its luminosity, could be underestimated by 
more than two orders of magnitude. If such massive black holes reside in galactic nuclei at the 
present epoch, they may be distinguishable from lower mass black holes accreting at rh > m c (e.g. 
Seyfert galaxies), primarily due to their hard X-rays and synchrotron radio emission. 

The proposed model could be roughly compatible with the observed X-ray luminosity 
evolution if (i) the transition in accretion flow occurs at L = 10 -2 Lem or to = m c = 0.3a 2 ~ 10~ 2 
or a ~ 0.3, (ii) ri n = 0.1 before transition and L cx to 2 after transition, and (hi) QSO remnants 
have black holes with M > 1O 9 M . These requirements are based on our particular set of initial 
conditions. For a more realistic model, some scatter in z% and/or 5 is inevitable, for which we 
lack detailed information. Despite these serious uncertainties, the simple model demonstrates 
a natural evolutionary trend which essentially needs only the second point (independent of our 
specific proposal) both in optical and X-ray. 

If QSOs at z ~ 2 — 3 are powered by black holes accreting at near Eddington rates, our 
explanation indicates m c ~ a few xl0~ 2 ora~ O(0. 1) where the latter value is interestingly close 
to those often quoted in Galactic accretion systems (Yi 1996, Narayan et al. 1996). The critical 
accretion rate in our model is related to several theoretical problems including the necessity for 
the two temperature plasma (Narayan & Yi 1995a, Rees et al. 1982, Phinney 1981), which merits 
further investigations. 

Alternatively, QSOs could be driven by electromagnetic power extracted from rapidly spinning 
black holes (Blandford & Znajek 1978). Although the spin-up of massive black holes from their 
formation is not clearly understood (cf. Thorne 1974), the spin-down of black holes does provide 
a time scale comparable to the QSOs' cosmological evolution time scale (Park & Vishniac 1988) 
when M ~ Mem- Such an evolutionary possibility still lacks an explanation as to why the 
accretion luminosity which is comparable to the electromagnetic power (Park & Vishniac 1994) 
also decreases on a similar time scale. 

The author thanks Masataka Fukugita, Ramesh Narayan, Martin Rees, Ed Turner, and Ethan 
Vishniac for useful discussions and comments. The author acknowledges support from SUAM 
Foundation. 
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Fig. 1. — Luminosity evolution index K(z) in the two cosmological models with a constant mass 
accretion rate for two different initial epochs {%i = 3,4). The curves, marked by the relevant 
parameters, describe the luminosity evolution for z < Zi and show effects of cosmological model 
and Zi. Evolutions at z < Z{ show a wide range of K including K ~ 2 — 3 with z-dependence. The 
slope is steeper for the flat universe (f^o = 1)- K(z) becomes invalid as z — > Zi as the accretion flow 
is expected to enter the high efficiency regime (see text). 
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Fig. 2. — Scale-free luminosity evolution of a QSO in the four cosmological models described in the 
text. The cosmological model corresponding to each curve is marked by the number and the critical 
redshift z c is shown in the figure. In this example, a = 0.3, (3 = 0.5, and Z{ = 3.5. For comparison, 
L(z) oc (1 + z) 3 is shown (solid line). The evolution becomes steeper than (1 + z) 3 at z < 0.5 as 
the mass accretion rate starts to reflect the exponential factor on a time scale t evo i = f3t age with 
= 0.5. 



